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Abstract 

Eight production quality PbW04 crystals of CMS have been irradiated in a 20 GeV/c proton beam up 
to fluences of 5.4 x 10^'^ cm~^. The damage recovery in these crystals has been followed for over a 
year Comparative irradiations with ^''Co photons have been performed on seven other crystals using 
a dose rate of 1 kGy/h. In proton irradiated crystals the light transmission band-edge shifts and the 
induced absorption length /iiND(A) oc A^^. In 7-irradiated crystals the band-edge does not shift but 
formation of absorption bands is seen clearly. The absorption length induced by 7-radiation in crystals 
verified to have excellent radiation hardness, saturates at a level below 0.5 m^ ^ . In the case of protons, 
we observe no correlation with the pre-characterised radiation hardness of the crystals and the induced 
absorption increases linearly with fluence. After a fluence of 5 x lO^'^cm^^, an induced absorption 
length of '^ISm^^ is seen with no sign of saturation. These observations provide strong evidence 
that high-energy protons create damage that cannot be reproduced with 7-irradiation. However, these 
hadronic effects manifest themselves only at integral fluences well beyond 10^^ cm^^ and most likely 
would escape undetected at lower fluences. A large fraction of the damage, both in proton- and 7- 
irradiated crystals, is either stable or recovers very slowly. 
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1 Introduction 



Lead Tungstate (PbWO 4) crystals will be used by several high-energy physics experiments QEIEI because they 
provide a compact homogeneous calorimeter with fast scintillation. In view of their use in the CMS experiment, 
where the calorimeters will face the harsh LHC conditions, the radiation hardness of PbWO 4 crystals has been 
extensively studied over the last decadej4] |5] IH [H |9l. These studies have mostly been done at 7-irradiation 
facilities, complemented with some neutron tests at reactors f lOV At LHC, however, calorimeters will also be 
exposed to a large fluence of high-energy hadrons, resulting from up to 8 x 10* s^^ minimum bias pp interactions 
at Vi=14TeV. 

Although the radiation damage caused by minimum ionising particles in PbWO 4 has been thoroughly studied 
by the 7-irradiations, this cannot be considered sufficient for LHC applications. In particular, the fundamental 
difference between few MeV photons and energetic hadrons - the effects of nuclear interactions in the crystal - 
has not been investigated in detail so far Some pion irradiations have been performed 1111 [T2I but no systematic 
study has been extended up to the full integrated hadron fluences expected at the LHC, although the neccessity of 
this has been advocated long ago 1 13J. 

Inelastic nuclear interactions break up the target nucleus and thus create impurities and distortions in the crystal 
lattice. These effects, however, are expected to be negligible at LHC fluences, since the natural impurity concen- 
tration of the crystals is relatively high. Another unique feature of hadronic interactions is the very dense ionisation 
of the created heavy nuclear fragments, which can have a range of up to 10 /im. Along their path they displace a 
large number of lattice atoms, but also ionise much more densely than a minimum ionising particle. In PbW04 
the fragments with highest ionising dE/dx are known to come from fission of lead and tungsten where the cross 
section for such reactions exceeds 100 mb 1 14J. Typical fission fragments, like Fe and Zr in Fig^ have energies 
up to 100 MeV. It can be seen from Fig.^that the dE/dx for such fragments is four orders of magnitude larger 
than for minimum-ionising protons - indicated by the solid dot in FigQ Since the thresholds to induce fission are 
several hundred MeV, tests with reactor neutrons or low-energy protons are not suitable to probe this regime. 




Figure 1: Simulated ionising energy loss of different fragments in PbW04. Non-ionising energy loss is not in- 
cluded, but even for heavy fragments it is negligible above ^10 MeV. This is quantified for lead fragments, where 
the total dE/dx is indicated by the dot-dashed line. The dashed circle shows the typical dE / dx-values of fission 
fragments at the beginning of their track. 

An irradiation test of PbWO 4 with hadrons is more complicated than a standard 7-irradiation. The production of 
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Table 1: List of all crystals, their initial radiation hardness pre-characterisation (/iJ'^*Q(420) in m~i) based on the 
standard CMS 7-irradiation procedure 1 16 J. ti„ is the exposure time in the proton beam with nominal flux (ftp or in 
the 1 kGy/h 7-facility, used in this work. 

nuclear fragments causes the crystals to become radioactive, up to a level which does not allow handling them for 
several weeks after exposure. Any front-end electronics or optical system, like photomultipliers, laser or fibre is 
also susceptible to damage in the hadron beam and would make the determination of an effect in the crystal itself 
more complicated. Therefore our aim was to make the simplest possible test we could imagine. We irradiated 
bare crystals and measured the light transmission after irradiation. Our study was targeted at stable or very slowly 
recovering damage after a proton fluence of up to 5 x 10i3 cm^^ 



2 The crystals 

For our studies we used PbWO 4 crystals produced by the Bogoroditsk Techno-Chemical Plant (BTCP) in Russia 
for the electromagnetic calorimeter (ECAL) of the CMS experiment 1 1 1. The crystals have the shape of truncated 
pyramids with nearly parallelepipedic dimensions of 2.4x2.4cm2 transversely and 23 cm length. All crystals were 
of production quality, and only because of slightly non-compliant mechanical dimensions or mild mechanical 
surface damage they could not be used in the ECAL construction. In particular, they were all pre-characterised 
for their radiation hardness through 7-irradiation at the Geneva Hospital *^°Co facility following a standard CMS 
qualification procedure. Throughout the paper we will refer to the results of this pre-characterisation as /^Jj'd (420) . 

The studies presented here concentrated on an examination of the change in longitudinal transmission (LT) at a 
wavelength of 420 nm, which corresponds to the peak emission of PbWO 4 scintillation light. The CMS technical 
specifications 1151 require the radiation-induced absorption coefficient 

at A=420nm to be below 1.5 m~i for a saturated damage from 7-irradiation at wlOOGy/h. In Eq.niLTg {LT) is 
the longitudinal transmission value measured before (after) irradiation and I is the crystal length (23 cm) through 
which transmission is measured. 

For the proton irradiations we had 12 crystals at our disposal. In order to minimise the effect of crystal-to-crystal 
variations, we selected for our irradiation 8 crystals which showed the most consistent behavior in optical charac- 
teristics and in terms of ii'l^'^{A2Q) . The selected crystals are listed as a-h in Tabled 

For the 7-irradiations, performed at a dose rate corresponding to IQi^ p/cm^, we used further 6 crystals, from a 
later production, which are listed as u-z in Tabled The pre-characterisation of these 6 crystals was done for 1 h 
at 350Gy/h which increases /i^Jp(420) by about 10% 1 16 1 with respect to 2 h at 250Gy/h, as used for the other 
crystals. Crystals u, v and w were of very good quality while the three others showed after the pre-characterisation 
quite high iiJ^'^{A2Q) values. Two of them (x and z) were actually outside of the CMS specifications. In addition 
to these 6, one of those remaining from the earlier 12, was considered good enough to be used in the 7-irradiations 
and is labelled as t in Tabled 

The target dose rate in the 7-irradiations was 1-2 kGy/h but an exact analysis of the actual rate is presented later 
Throughout this paper we use the convention to assign small letters to crystals irradiated at a dose rate of ~ 1 kGy/h 
and capital letters to those irradiated at rates of 10 kGy/h. A prime or a double-prime after the crystal identifier 
indicates a second or third irradiation of the same crystal. 
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The <iii times given in Table[2for the proton irradiations are the target values but the actual fluences were deter- 
mined by standard activation techniques. For the 7-irradiations, performed in static and very stable conditions, the 
iirr values are the actual exposure times which are used in the dose determination. 

3 Irradiations 

For the proton irradiations we used the IRRADl facility 1171 in the T7 beam line of the CERN PS accelerator, 
which at the time of our first irradiations was delivering a proton beam of 20 GeV/c and during a later campaign 
24GeV/c. In order to obtain a uniform irradiation over the whole crystal front face, we requested a special beam, 
which had not been produced before in that facility. In particular, we asked for the lowest beam intensity which 
the PS could provide in a controllable way in the T7 beam Une. 

At 20GeV/c we obtained a flux of ~ 10^^ cm^^h^^ over an area of roughly 4.5 x 4.5cm^. At 24GeV/c a further 
intensity reduction by a factor of 2 was reached. The uniformity was not perfect over the whole beam spot, but a 
fairly uniform area covering the crystal could be found. In this low-intensity mode the PS was delivering on average 
one proton spill every 42 seconds. The duration of the spill was 600 ms at 20 GeV/c and 400 ms at 24 GeV/c. 

In order to investigate a possible rate dependence, we performed some irradiations at a higher intensity, of ~ 
lO^'^ cm^'^h^^. To achieve this, the spill frequency was increased to 2 spills every 14s and the beam spot was 
reduced to a size of about 3x3 cm^ which was still sufficient to cover the crystal while maintaining an acceptable 
beam uniformity. 

The crystals were brought into the beam by a remotely controlled shuttle. The beam was hitting the small front 
face of the crystal and was parallel to the long crystal axis. 

The 7-irradiations were performed at the Calliope plant of ENEA-Casaccia, Italy. The facility provides a ^°Co 
source which at the time of our irradiations had a total activity of 750 TBq. The 48 source rods were arranged 
in two concentric cylinders of about 20 cm outer radius and 26 cm height. We irradiated simultaneously up to 5 
crystals, standing side-by-side at a distance of 40 cm from the source centre, with their long side parallel to the 
source axis. Each crystal was housed in its own styropor box. The distance between two adjacent crystals was 
about 34 mm. In order to obtain as uniform an exposure as possible, each box was turned once by 180° around its 
long axis at half-time through its exposure. Unfortunately the crystals could not be precisely adjusted with respect 
to the source and an offset of 52 mm remained along the axis. Thus one extremity of the crystals was slightly less 
exposed. 

4 Measurement methods 
4.1 Light transmission 

The LT measurements were taken with a Perkin Elmer Lambda 900 spectrophotometer for wavelengths between 
300 and 800 nm in 1 nm steps. Since the grating which allows to select the desired wavelengths polarises light, we 
used a wide-band depolariser to be reasonably free of polarisation effects. The width of the beam is defined by 
the slit opening which fixes the wavelength distribution: this we set to AA = 1 nm which yields about 7 mm beam 
size. The height was adjusted with a "common beam mask" to approximately 10 mm. 

We will discuss later the implications of non-uniform proton beam profiles, which required us to align the crystals 
in a reproducible way. For each measurement we checked that the entry and exit points of the light were centred 
on the crystal face. The accuracy of this alignment was roughly ± 1 mm. Measurements with intentionally larger 
displacements showed this to be sufficiently good not to introduce significant variation for any of our crystals. 

For each proton-irradiated crystal, at least one set of LT measurements was performed, where the light was passed 
not only through the centre, but also as close to each of the four sides as possible. The information gained from 
these "scans" served the purpose of understanding effects of beam non-uniformities. 

Before each session of measurements we experimentally verified the absolute calibration of the device by measur- 
ing the LT of a non-irradiated reference crystal. Fig.|2] shows the results of these measurements as a function of 
time at 420 nm. From the scatter of the points we derive a la deviation of 0.58%. However, close to the band-edge 
the transmission changes rapidly and the accuracy of the wavelength dominates. We were able to reproduce the 
scatter of the reference crystal LT close to the band edge, if we assumed that the wavelength is reproduced with 
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Figure 2: LT data measured on tlie reference crystal at 420 nm. The 0-point of the time axis is arbitrary, but the 
relative time intervals correspond to real elapsed time. 

0.23 nm accuracy. Thus we estimate the total error of the LT measurement at wavelength A as 

A[LT(A)] = 0.0058 x LT{X) ® 0.23 x 6lt{\}, (2) 

where ^lt(a) is the difference between two adjacent LT points, 1 nm apart. The second term is significant only 
close to the band edge, when Si^t{X) can be as big as three percent-units per nm. At large wavelengths, when the 
transmission changes slowly with A, the error is given by the first term alone. 



4.2 Induced radioactivity 

The induced radioactivity dose rate (i^ind) was measured with an Automess 6150AD6 1181 at a distance of 4.5 cm 
from the long face of the crystal at its longitudinal centre. The device was calibrated to yield the photon equivalent 
dose, which is the quantity extracted from the simulations described in appendixlAl 

The reference point of the sensitive element in the 6150AD6 is reported to be 12 mm behind the entrance window. 
Thus our actual distance was 57 mm from the crystal face. The -Dind values of the least active crystals were 
comparable to natural background, which we determined by integrating over a long period with the 6150AD6. In 
order to verify the linearity of the device even at low counting rates, we also measured the dose from a standard 
point-like ^■^^Cs source at various distances. Both methods gave a background of 0.06 /iSv/h and we observed no 
dependence on counting rate. 

The nominal energy range of the 6150AD6 is 60keV-L3MeV. Below 60keV the efficiency drops abruptly, but in 
our crystals only a negligible fraction of dose comes from low energy photons. The efficiency curve 1 18 1, which 
rises steeply above 1 MeV, is cut at 1.3 Me V, which marks the end of the calibrated range. Our assumption of 
constant efficiency above this energy most likely leads to an underestimate of our simulations with respect to the 
measured values. 

Since PbW04 attenuates photons quite efficiently, a measurement of Diud provides a means to approximately 
determine the proton beam uniformity. For this purpose Diud was measured on all long sides for each crystal. 
Although the attenuation is exponential, a simple average over the four sides gives the total activation with an 
accuracy of a few percent even for our most non-uniform irradiations. 

All Dind measurements were prolonged sufficiently to reduce the statistical uncertainty to the level of 1-2%. 



5 Simulation of the irradiation conditions 
5.1 Proton beam 

Since the crystals are longer than one hadronic interaction length, the protons initiate an intense hadronic cascade. 
In addition, the irradiation zone is a fairly small volume with steel walls on 4 sides at distance of few tens of 
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Figure 3: Ionising dose, star density, charged hadron and neutron fluences per unit of incident proton fluence 
as obtained from the simulations. Values are for a 5x5cm^ proton beam of 20GeV/c (solid dots and lines) or 
24GeV/c (open dots, dashed lines). The lines and associated values give the average over the whole crystal length 
in the central area, covered by the light beam in the LT measurements. 
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Table 2: Simulated hadron fluxes (cm ^ per unit of incident flux) in the first and last millimetre of the 23 cm long 
crystal hit by a 20 GeV/c proton beam. 

centimetres. The back wall is in fact the beam dump of the T7 beam line. Thus the possible effects of back- 
scattered radiation cannot be ignored. 

In order to properly understand the radiation field within and around the crystal - and in particular the influence 
on the determination of the incident proton fluence - a full simulation of the irradiation setup was performed with 
the FLUKA code 1 19 1. For the simulations, the crystal shape was described as a 23 cm long truncated equal-sided 
pyramid with a minor face of 22 x 22mm^ and a major face of 26x26mm^. The main quantities characterising 
the radiation field, as obtained from these simulations, are shown in Fig.|3l 

Damage caused by 7-irradiation is found to saturate at a level corresponding to the dose rate used II20I . How- 
ever, since we are looking for the possible existence of cumulative damage^^ specific to hadrons, it seems more 

^' Earlier hadron irradiations reported in IlllfT^ l did not show any difference between 7 and pion irradiation performed at 
the same rate. But these tests cannot be considered sufficient to exclude possible cumulative effects because they were not 
extended to high enough fluences. 
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Figure 4: Energy spectra of protons, charged pions and neutrons in the first (soUd circles) and last (open circles) 
millimetre of a 23 cm long crystal. The dotted line shows the 20MeV threshold above which hadrons can produce 
stars. 

appropriate to use, instead of dose rate, a purely hadronic quantity against which we could plot the damage. 

The incident proton beam is practically monoenergetic but, due to cascade formation, the hadron spectrum within 
the crystal, shown in Fig.@] comprises all energies from sub-MeV neutrons to full beam energy. Table|2lshows the 
integral fluences at beam energy and above 20MeV - the threshold indicated by the dashed lines in Fig.|4] With 
such a wide energy spectrum and since a dependence between damage and hadron energy cannot be excluded a 
priori, it is not obvious how damage could be related to hadron fluence. 

It is plausible that any damage unique to hadrons should be caused by the component of the hadron spectrum which 
is producing inelastic interactions. The most suitable quantity to describe the radiation field for our purpose, is 
called star density, which was originally introduced to parametrise radioactivity induced by hadron irradiation 1 2 1 1 . 
A star is defined as an inelastic hadronic interaction caused by a projectile above a given threshold energy. This 
means that the star density is actually the integral of the total hadron fluence above a given threshold, weighted by 
the inelastic cross section. The star density is exactly zero for irradiations with low energy photons or neutrons 
below the threshold energy, which we take as 20 Me V. For constant beam conditions the star density is of course 
proportional to the fluence. In our analysis we will prefer the latter because it can be experimentally determined 
with good accuracy. The simulated star densities, however, will be used to scale between different beam energies 
and to compare different radiation environments. 

Since the radioactivity produced in the crystal is roughly proportional to the rate at which stars are produced 1221 . 
a measurement of -Dind in a crystal can provide an estimate of the star density. This is discussed in more detail in 
appendixIXl 

5.2 ^°Co source 

The dose rates at the Calliope plant 1123 1 are well known from measurements, performed regularly by the operators. 
However, these values refer to dose in air or water in the absence of scattering material. In order to determine 
the dose profile within the crystals, detailed simulations of the irradiation setup were mandatory. Fig.|5lshows the 
FLUKA simulation model of the source and 5 crystals in their boxes. Since durations of irradiations were different, 
not all positions were occupied at all times. Due to these changes in configuration and the slightly different position 
with respect to the source, the dose in individual crystals is slightly different. 

As in the case of proton irradiated crystals, we again determine the transmission along the long axis of the crystal, 
and therefore we are interested in the dose along the centre. Fig.|6lshows the simulated transverse dose rate profiles, 
which exhibit a pronounced U-shape, caused by turning the crystal half-way through the irradiation. The fairly flat 
central part, of about ±5 mm width, fits very well the width of the light beam in the LT-measurements. In the other 
dimension the profile is flat. 
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Figure 5: Top view of the simulation geometry representing the Calliope facility. The two concentric assemblies 
of source rod rings are seen on the right and the 5 crystals inside their boxes on the left. The positions of each 
crystal are indicated by their identifiers next to the boxes. The small dots on and within the crystal boxes represent 
alanine dosimeter positions. 

Figure0shows the simulated longitudinal dose rate profile along the long axis of the crystal. The dose drops quite 
significantly towards the end of the crystal because the height of the source could not be exactly aligned with the 
crystals. The discontinuities, e.g. the step at 18 cm for crystal w, are caused by the support structures of the source 
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Figure 6: Simulated transverse dose rate profiles in the crystals. Shown are the minimum (crystal t) and maximum 
(crystal w) dose rates in the 7 crystals. The profiles represent an average over the whole length of the crystal. 
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Figure 7: Simulated longitudinal dose rate profiles along the centre of the crystals (±5 mm around central axis). 
Shown are the minimum (crystal f) and maximum (crystal w) dose rates in the 7 crystals. 

assembly. 



6 Fluence and dose determination 
6.1 Proton fluence determination in IRRADl 

The fluence, <i>p, of the incident proton beam was determined with aluminium foils, cut to a size of 24 x 24 mm^ 
and placed in front of the crystals at a distance of about 10 mm. Secondaries produced by the crystals, neutrons 
in particular, influence this measurement. Therefore we preferred to use the isotope ^^Na instead of the more 
standard ^"'Na. The determined fluence, in particular its relation to the star density within the crystal, is sensitive 
to the geometry and alignment. A further factor to be taken into account came from a change of beam momentum 
from the initial 20 GeV/c to 24 GeV/c. Therefore we paid great attention to a careful fluence analysis and cross 
checks with the induced radioactivity in the crystal itself. The technical aspects of this and all effects considered, 
together with the correction factors derived by us, are detailed in appendixiBl 

Table|3l summarises the proton fluences and gives the dose rates and total doses averaged along the axis of the 
crystal. The first column, <i>p(^^Na), gives the fluence value obtained from the aluminium foils but corrected by 
the factor of 0.89 in order to account for activation by secondaries. The ^ (l>p-column gives the integral of all 
subsequent irradiations for a given crystal and also includes the correction by a factor of 1 .086 to account for the 
higher beam energy in the case of crystals a" and h. Any quantity that is cumulative and stable, should appear 
proportional to ^ $p. 

The dose values given in Table|3lcorrespond to an average along the axis of the crystal, i.e. along the path of light 
in the LT measurements. Averaging was done over the dimensions of the light beam. It can be seen from Fig.|3l 
that in the proton irradiations the maximum dose rate, at a depth of about 10 cm within the crystal, is ^^30% above 
this average. 

One issue which will be discussed in detail in appendixicl is the uniformity of the beam spot across the crystal 
face. These effects, however, remained small enough to require no corrections. Actually, as shown in appendixicl 
the simulated effects of non-uniform irradiation agree nicely with measured variations of Di„d and LT close to the 
edges of the crystal. 
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Table 3: Proton fluences $p(^^Na) for each individual irradiation as obtained from the Al-activation and cumulative 
fluences (^ (f>p) for a given crystal. The latter also includes the correction for beam energy (see text). The last 
two columns give the average ionising dose rate and integral dose along the axis of the crystal. The dose values 
are based on simulations, renormalised to the measured fluence at proper beam energy. 

6.2 Dose determination at the 7-facility 

Alanine dosimeters provide a well-established method to determine cumulative ionising doses. The absolute cali- 
bration is usually done with a standard ^°Co source and converted to give the dose in air. Ten alanine dosimeters 
which accompanied the irradiations provided a crucial cross-check of our dose simulations. Six of these dosimeters 
were fixed on the "front" of each of the crystals u-z, i.e. the side facing the source, while two were on the "back" 
of crystals v and y. When turning the crystals, "front" and "back" got inverted so that each of these eight alanines 
was averaging equally over the "front" and "back" of a crystal. A further two dosimeters were fixed on the source 
side of the styropor boxes housing crystals v and y. These dosimeters were removed when the boxes were turned 
after the first 23 hours of exposure. All 10 alanines were included in the simulations and the dose in them was 
recorded. Figure|8lshows a comparison between the doses determined from the alanines and the simulated values. 
While the low doses of the 1 h and 8.3 h irradiations are in prefect agreement, alanines exposed for the longest 
times indicate a lower value than the simulations. Due to the symmetry of the simulated setup - essentially the 
same result was just rescaled with time - it is inconceivable that the simulation would fail for higher doses. The 
response of alanine is usually quoted to be linear up to about 100 kGy. However, when approaching the upper 
limit, saturation sets in gradually. Although saturation is taken into account in the calibration 1241 . an explanation 
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Figure 8: Simulated (open circles) and measured (solid circles) cumulative doses in the 10 alanine dosimeters 
exposed during our 7-irradiations. 
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Dose (kGy) 


u 


1.02 


1.02 


z 


1.15 


1.15 


X 


1 09 
i .uz 




w 


1.15 


9.59 


t 


1.00 


27.3 


y 


1.06 


48.9 


V 


1.09 


50.3 



Table 4: Dose rates and cumulative doses reached in the 7-irradiations. Values are averages over the crystal length 
along the centre axis. The statistical errors are smaller than 1%. 

for the difference observed could still be if saturation depends also on the photon energy spectrum, which is softer 
in the presence of the crystal than for a plain ^°Co calibration source. 

Table|3 summarises the doses and dose rates reached in the 7-irradiations. These values have been extracted 
directly from the simulation described in section]^] but their accuracy has been experimentally verified by the 
comparison with the alanine data, shown in Fig.ls] The dose rate at one end was about 10% higher than the average 
and dropped to 20% below average at the other end of the crystals, as can be seen in Fig.0 



7 Results 

Fig.|9] shows typical LT curves of a non-irradiated crystal and of some crystals which have suffered radiation 
damage in the proton beam. With increasing ^ $p, not only the transmission worsens, especially in the blue 
region of the spectrum, but in addition the LT band edge shifts towards higher wavelengths. At ^ <l>p « 10^^ cm^^ 
this shift starts to increasingly cut into the intrinsic emission spectrum of the crystal. 

Fig-CSlshows the transmission curves as a function of wavelength for all 7 crystals irradiated with ^''Co photons. 
The shapes of the LT curves of 7-irradiated crystals are completely different from those of proton-irradiated ones. 
Even after the highest cumulative dose reached, the band-edge in 7-irradiated crystals does not shift at all, thus 




Wavelength (nm) 



Figure 9: LT curves for crystals with various degrees of proton induced radiation damage. The emission spectrum 
(dotted line) is taken from 1 25 1 and has arbitrary normalisation. 
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Figure 10: LT curves for crystals with various degrees of 7-induced radiation damage. The lines show the light 
transmission prior to irradiation. The order of these curves is the same as of the post-irradiation ones. The emission 
spectrum (dotted line) is taken from 1 25 1 and has arbitrary normalisation. 
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Figure 11: The fits of Eq.|3lto recovery data of our proton irradiated crystals. The last points of crystal G are 
excluded from the fit - see text. 

establishing a fundamental difference between the damage caused by protons and the one caused by purely ionising 
radiation. Instead, one recognizes the absorption band around 420 nm, which is typical for 7-irradiated BTCP 
crystals |5 1. Although this band is probably also present in proton-irradiated crystals, it does not appear in the LT 
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Figure 12: Recovery data of crystals irradiated with ^"Co photons up to various cumulative doses at a dose rate of 
-l.lkGy/h. 

curves as a main feature of the damage. Furthermore, in 7-induced damage we see a significant crystal-to-crystal 
variation. This agrees well with the variation of the /i^^*^(420) values in Tabled No significant correlation of this 
kind is seen in the proton-irradiated crystals. 



7.1 Recovery time constants 

We fit the time dependence of /xind with a sum of two exponentials and a constant 



Mind (420, i 



/ t 



(3) 



where t^cc is the time elapsed since the irradiation. Attempts to fit Eq.|3]for each crystal, with all 5 parameters free, 
revealed that the longer time constant (T2) was comparable to, or larger than, iroc and thus the correlation between 
the second term and the constant was too strong to allow for a reliable fit. It is, however, reasonable to assume that 
the damage corresponds to the same color centres, each of which corresponds to a unique time constant. Under 
this assumption we fitted the for all crystals simultaneously and, with these fixed, the amplitudes Ai for each 
crystal separately. This process was iterative, i.e. Tj and Ai were fitted repeatedly with increasingly better initial 
guesses until convergence was reached. The quality of the fits, in terms of x^, was comparable to that obtained by 
the 5-parameter fit of Eq.l^for each crystal separately. 

Figs.^2and[2]show the fits for proton and 7-irradiated crystals, respectively. In comparison to the observed 
recovery time constants, crystals a, a ' and F were re-irradiated too soon and therefore do not allow to obtain 
reasonable fits. 

Of the proton irradiated crystals, d was annealed 240 days after irradiation, thus we have no longer-term room 
temperature recovery data for it. While essentially no recovery of crystal d was observed at 160 °C, partial recovery 
took place at 250 °C. This temperature was sufficient to restore the pre-irradiation transmission at A >400 nm. At 
a temperature of 350 °C almost complete recovery was reached, except for a small residual shift of the band-edge. 

In the case of crystal G a peculiar step appeared after 250 days, which we traced back to an unintentional exposure 
to normal room light from fluorescent tubes for about 60 hours during a long radioactivity measurement. The 
existence of such light-induced annealing was qualitatively verified by exposing to the same room lights a crystal 
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Figure 13: The amplitudes Ai of Eq.|3]as a function of dose. The dose to fluence correspondence for the proton- 
irradiated crystals can be found in Table|3l When symbols are not present in the plot, the fitted amplitude for that 
crystal is consistent with zero. 
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Figure 14: Plots of Eq.|3lwith the parameter values shown in Fig.^] after renormalisation to unity at t=60days. 
The symbols indicate the range over which data were available for a given crystal. 

that had been proton irradiated while exercising the procedures for this work. The very severe damage in that 
crystal had remained stable for almost 2 years but about 15% of it annealed fairly rapidly under light exposure. 

The 7-irradiated crystal w was annealed 65 days after irradiation, in preparation for a later proton irradiation. 

From the global fit we obtained Ti=17.2days and T2=650days. The amplitudes corresponding to these recovery 
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time constants are shown in Fig.[0] The amplitude corresponding to ri is slightly increasing with dose in the 
case of proton irradiation while it decreases as a function of 7-dose. This decrease could be due to the longer- 
lived defects, which build up and contribute to light attenuation, making Ai relatively less effective. Since T2 was 
fitted to a value of 650 days, which is longer than the period of the recovery follow-up, the amplitudes A2 and ^3 
could not be reliably resolved from each other This was particularly true for all 7-irradiated crystals and proton 
UTadiated crystals d and G, where recovery data extends only to tree ^200 days. However, the sum of A2 + A3 for 
the proton irradiated crystals increases linearly with fluence, while it is consistent with a constant for 7-irradiated 
crystals. 

All fits - re-normalised to the same values at ii.cc=60days - are shown in Fig.^^ Re-normalisation at 60 days 
was chosen such that the ti -component of damage has practically disappeared. Thus the left side of the re- 
normaUsation point shows the relative amplitude of the damage corresponding to ri, while the right side provides 
a visual comparison of the magnitude of the long-term damage. Especially in terms of the short component the 
crystals group according to the integral fluence or dose they received. Those exposed to J2^p ^ 10^^ cm^^ (or 
~1 kGy) recover, in relative terms, more than those irradiated to higher fluence. While a slight difference between 
proton- and 7-irradiated crystals can be seen in the magnitude of the fast component, the long term recovery 
appears to be independent of the radiation type. 



7.2 Correlation between light transmission and ^ $p or 7-dose 




Figure 15: Induced absorption /iind(420) as a function of cumulative proton fluence. 

FigurefTsl shows ^ind(420) as a function of fluence 150 days after irradiation as calculated from Eq.|3l using the 
parameter values shown in FigfTsl 

For crystals irradiated several times, the cumulative fluence ^ $p has been used and /iind (420) has been nor- 
malised to the transmission before the first irradiation. We have fitted the data, which cover two orders of mag- 
nitude, on a log-log scale and obtain an exponent that corresponds to only 1% deviation from linearity per order 
of magnitude in fluence. In particular, if the data in Fig. [21 are compared with Tabled no correlation between 
proton-induced damage and /i^ND (420) can be observed. 

Fig. II 61 confirms that even at a dose rate as high as 1 kGy/h the 7-damage saturates. After 1 hour at 1 kGy/h, 
Mind (420) has not yet reached its plateau value, but between 8.3 and 46 hours the differences observed can be 
entirely attributed to crystal-to-crystal variation. This becomes particularly clear when we normalise ^ind(420). 



14 




S a 

1.8 

^1 1.6 
^ 1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




I 



■ 


Crystal u 




Crystal t 


□ 


Crystal z 


▲ 


Crystal v 


• 


Crystal w 


O 


Crystal y 


A 


Crystal x 







o 
▲ 



10 15 20 25 30 35 40 45 50 

t. . (h at 1 kGy/h) 



Figure 16: Values of /iind (420) measured on the 7 7-irradiated crystals after 20 days recovery at room temperature. 
The upper plot shows the measured data and the lower one the ratio to the /ij^'o 

(420)-values from Table[T] 

as measured by us, with the /xJ'^'q (420) -values from Tabled which is supposed to cancel any crystal-to-crystal 
variation. It should be noted, however, that the iij^^{420) values are obtained only about 1 hour after irradiation, 
i.e. before short-term recovery is completed. In fact, it has been observed that about 60% of the damage con- 
tributing to /i^^'Q(420) recovers during the first 20 days while the rest is much more stable |16|. If this is taken 
into account in the /xind (4:20) /h]^^ (420) ratio of Fig.^] the 1 h points are nearly consistent with unity while the 
ratios for longer irradiations get pushed to a value of about 3. This large ratio is at least partly - maybe even mostly 
- due to the fact that the pre-characterisation irradiations did not last long enough to saturate all color centres. 

Crystal f shows slightly higher /iiND(420) values than the other 6 crystals in Fig.^] The deviation is even more 
pronounced in the /xind(420)//x^J'^(420) ratio. These differences might well be, because t is from the earlier set 
of 12 crystals received for the proton tests. These had undergone a different pre-characterisation procedure almost 
a year earlier. While the 10% difference in /ij^^*p(420), due to the different dose rate, has been taken into account, 
there could still be other systematic or procedural effects that might have led to a slightly lower /i/^o (420) of 
crystal t. 

When comparing the data of Figs. and ^] it should be taken into account that the proton irradiated crystals 
all had /xj'^'p(420) < 1 m^^ and most of them should be compared with the solid symbols in the upper plot of 
Fig.[^ In particular, as can be seen from Tabled crystal a" was characterised as having the same, excellent, 
radiation hardness as crystal v. After exposure to comparable integral doses, the proton irradiated crystal a " shows 
an induced absorption of 15 m^^ while the 7-irradiated crystal v has saturated at a value of only 0.3 m^^. 

The shape of the LT curves, in particular the shift of the band-edge after proton exposure and the fact that two 
universal time constants fit the data well but give very different dose-dependence for the amplitudes in Eq.|3 leads 
us to the following interpretation: 
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Figure 17: Plot of /iiND(A) against A for proton and 7-damaged crystals exposed to different fluences and doses. 
The dot-dashed line shows A^**, fitted to the proton damage curve. The good agreement is an indication of Rayleigh 
scattering from small centres of severe damage. 

• The damage corresponding to the two values is produced by proton and 7-irradiation in the same way and 
results in the same recovery-characteristics. Most likely this damage corresponds to the same color centres. 

• The band-edge shift is unique to proton irradiation and corresponds to a different damage mechanism. This 
damage appears to be stable and cumulative, at least over the ^ $p range explored. 

A similar band-edge shift has been observed earlier in proton-irradiated BGO-crystals flE] while no such shift was 
created by 7-exposure. A later re-analysis of the data in Ii26i revealed that the proton induced damage is cumulative 
also in BGOl271. 

The amplitudes Ai and A2 exhibit some dependence on proton fluence, which is not perfectly consistent with the 
above assumption. It is conceivable, however, that other defect centres, specific to protons, play some role and 
have time constants that could not be resolved from our two universal values and the constant term. 

The assumption that the highly ionising fragments, produced by inelastic hadronic interactions, would be respon- 
sible for the observed specific proton damage implies that there might be very small regions of severe damage 
distributed within the crystal. Fig.ll7l shows /iiND(A), plotted as a function of the wavelength. Above '^420nm 
the dependency for heavily proton-damaged crystals goes like A^'', which is exactly what would result if the at- 
tenuation was entirely due to Rayleigh-scattering. In 7-damaged crystals no such hnear dependence appears and 
only a structure due to absorption band formation is seen. The significant crystal-to-crystal variation of 7-damage 
is again evident. The Ijl'^^^{420) values in Table[2 suggest that the proton-damaged crystals should rather be 
compared with the less damaged 7-irradiated crystals u, w and v. For proton damaged crystals, the A^^ slope of 
Rayleigh-scattering is clearly superimposed on top of the absorption-band structure from pure ionisation damage. 
This provides firm evidence that the specific proton damage is of a very different nature than what can be caused by 
purely ionising radiation. At fluences of the order of 10^^ cm^'^ and below, the A"'' slope starts to get submerged 
in the absorption band structure which means that at low fluences it would be difficult to distinguish the specific 
proton damage. 

The /iiND (A) slope starts to deviate strongly from the A^** behaviour at wavelengths below 420 nm where the crys- 
tals become strongly absorptive. The Rayleigh scattering formalism is known to break down when the dimensions 
of the scattering centres become larger than ^10% of the wavelength. This means that the created defect centres 
must have dimensions ^0.4 fim. 

The fact that proton damage seems to be connected to very small, highly damaged, regions strongly points in the 
direction that the origin is indeed in the heavily ionising fragments created in inelastic hadronic interactions. In this 
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case, the overall dose rate of the irradiations should not be significant, since the local instantaneous ionisation of 
the fragment is overwhelmingly large. Although the PS beam did not allow to lower the intensity in order to obtain 
firm experimental proof for this, all evidence supports the assumption that the proton damage will be cumulative 
irrespective of the dose rate. 



8 Extrapolation to CMS at LHC 

The high-precision ECAL will play a crucial role in the CMS experiment for the search of the Higgs boson and 
exploration of possible new physics beyond the Standard Model. The calorimeter comprises ^76000 PbW04 
crystals, which provide hermetic coverage up to a pseudorapidity 77 = 3.0. In order to fully profit from the 
advantage of a homogeneous calorimeter, namely a very low stochastic term in the energy resolution, the constant 
term must be kept low as well. This can be achieved only by accurate calibration and permanent monitoring of the 
crystals. Radiation damage, not compensated by monitoring, would quickly degrade the calorimeter performance. 
Therefore, the LT of each crystal will be calibrated with a 440 nm laser. In addition Z e+e^ and W ^ ev 
events can be used for in-situ calibration. However, if the transparency of the crystals decreases, this will ultimately 
affect the photostatistics and thereby increase the stochastic term in an unrecoverable way. 

The endcaps (EE) are the parts of the ECAL which are exposed to the harshest radiation environment. The area 
where high resolution is desirable extends to |?;| = 2.5. The hadron spectrum incident on the EE varies with 
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Figure 18: Star density at different 77 in the CMS endcap ECAL for 500 fb ^ integrated luminosity. The histogram 
shows the average star density in a single crystal irradiated with 10^^ protons per cm'^ in the 20GeV/c beam-line. 

radius, but it is always broad and on average much less energetic than the 20GeV/c proton beam used in our tests. 
In Fig.^]we compare the simulated star densities at various ?;-values in the EE and along the centre of a crystal 
exposed to lO^'^p/cm^ at 20GeV/c. In terms of pure star density, a fluence of $p = 10^"^ cm^^ of 20GeV/c 
protons would correspond to the star density expected at 77 « 2.2 over 10 years LHC operation. However, the 
characteristics of stars change with energy. While in our beam test a majority of stars are generated by high- 
energy protons, in the EE they are predominantly due to pion and neutron interactions at energies rarely exceeding 
1 GeV. At low energy typically a few nucleons are emitted and one heavy slow recoil is produced. With increasing 
energy, the multiplicity of secondary particles increases and the mean atomic mass of the heavy nuclear fragment 
decreases, while its energy increases slowly. Energies of fission fragments are fairly constant but the cross section 
is energy-dependent. These effects can be taken into account only by a detailed simulation of the hadron-nucleus 
events followed by transport of the fragments down to zero energy. However, prior to such a simulation, it has to 
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Figure 19: Average sum of tracklength per star with ionisation density above a given value of (dE/dx)(.rit for the 
beam test and for CMS endcap ECAL conditions. 

be decided which is the physical process causing the damage. 

We have argued that stars are likely to be at the origin of the specific proton damage and our experimental results 
- especially the consistency with Rayleigh scattering - seem to support this. However, our irradiation data do not 
tell us by what mechanism the stars cause the observed damage. In high-resistivity silicon the hadron damage is 
known to be due to the created impurities and, especially, lattice defects 1281 . However, since the natural impurity 
concentrations in PbW04 are much higher, a similar damage mechanism cannot be assumed without a quantitative 
consideration. From Fig.|3lwe deduce an average of 0.14 stars/cm'^ per unit of flux. Thus, for $p — 10^'^ cm^^ 
we obtain an impurity concentration of about 3 x 10~^^ per atom. This is far lower than the natural density of 
imperfections in the crystal, which is at the 10~^ level. 

However, along their path, the recoils collide elastically with atoms of the lattice and kick these out from their 
position. Most of the displacements heal immediately, but some remain stable or metastable. Simulation methods 
of fragment transport have been developed for studies of Single Event Upsets in silicon based electronics chips 1 29 1 . 
Thereafter they have been successfully applied to bulk damage studies of silicon detectors 1 28 1. Silicon, in contrast 
to PbW04, is a well known material from this point of view. Thus the adaption of the simulation to the crystals 
introduced several uncertainties. In particular, the displacement energy, about 20 eV in silicon, is not known 
for PbW0 4. While it is likely to be higher, it is unlikely to be much different from the silicon value. For our 
simulations we assume a value of 50 eV. The simulations proceed as described in 1 28 , 29 1 and give on average 
2700 displaced atoms per 20 Ge V/c proton. Counting both the vacancies and the interstitials as defects, the number 
doubles. However, at least 90% of the vacancy-interstitial pairs will recombine immediately, so a reasonable 
assumption for the number of defects created per proton is <500. This gives an upper limit of ~ 10^^ for the 
defect concentration per lattice atom, which still remains about two orders of magnitude below the natural impurity 
levels. 

A more plausible damage mechanism seems to be that the created nuclear fragments cause extremely high local 
ionisation densities along the short path which they cover in the crystal. It can be seen from Fig.^that a lead 
ion can reach a maximum dE/dx of 5 x 10^ MeV/cm, which is 50000 times larger than for a minimum ionising 
particle. It is certain, however, that Pb fragments with 1 GeV kinetic energy will not be created by hadron-nucleus 
collisions. The most energetic fragments are due to high-energy induced fission of W and Pb and typically have 
energies of 100 MeV or just below. Zr and Fe, shown in Fig.Q are good examples of possible fission products. 
Thus the maximum dE/dx of the fragments can be about 2x10'' MeV/cm. If these very densely ionising fragments 
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are responsible for specific hadron damage, then the main question is what is the critical value {dE / dx)crit needed 
to create a stable damage centre. If such a value can be identified, then the amount of damage is quite likely 
to scale with the total tracklength of fragments exceeding (dE/dx) crit- This tracklength is shown in Fig.ll9lfor 
(dE /dx)ci-it values above 1 00 Me V/cm. The tracklengths are normalized per star, which allows a direct comparison 
of the reactions induced by 20 Ge V/c protons and those expected in the CMS EE, for which an average spectrum 
between 77=1.9-2.8 has been used in the simulation. In order to be insensitive to the non-uniformity close to the 
crystal edges, the tracklengths were recorded only ±5 mm around the centre line of the crystal, which corresponds 
to the region mapped in the LT measurements. 

The difference in tracklength between the 20 GeV/c proton induced events and those created by the EE spectrum 
shows only a weak dependence on (dE/dx)crit and amounts on average to a factor of 4. At {dE / dx) be- 
low 1 GeV/cm, the tracklength is mostly due to slow protons. For these it should be noted that the simulations 
consider only protons produced by the stars. Especially in the CMS EE, but to some extent also in our beam 
test, protons produced elsewhere can enter the crystal and produce a highly ionising track upon stopping^'. For 
(dE / dx) crit >1 GeV, heavier fragments dominate because, as shown in Fig.^ the dE/dx of a proton can never 
exceed 1.5 GeV/cm. It should be pointed out that the results of Fig.^]are not sensitive to our assumption of 50 eV 
as the displacement threshold. That parameter affects only the total number of produced lattice defects, but not the 
highly ionising tracks. 

With the provision that these simulations are based on an unproven assumption of the mechanism causing the 
specific hadron damage. Fig. ^] predicts that stars produced by our beam test are on average a factor of 4 more 
damaging than those expected in the EE. With 20 GeV/c protons the specific hadron damage becomes dominant at 
J2^p ~ 10^'^ cm^^, when /iiND(420) w 2m^^. Taking the factor of 4 into account but assuming that the damage 
is stable over time, this would correspond to roughly 77=2.6 after an integrated LHC luminosity of 500 fb~^. If this 
extrapolation is correct, the performance of the high-resolution region of the ECAL will not be compromised by 
hadron damage during the first 10 years of LHC operation. 

9 Conclusions 

We have performed an extensive irradiation campaign, aimed at verifying the possible existence of specific, hadron- 
induced damage in PbW04 crystals. Emphasis has been put on a careful analysis and detailed understanding of 
the irradiation conditions and a long follow-up of damage recovery. 

The very intense radiation environment would have damaged any front-end electronics or other sensitive equip- 
ment, making a determination of damage on the crystal itself complicated if an in situ determination of LT had been 
attempted. To avoid such ambiguities we put our emphasis on a simple procedure aimed solely at the determination 
of long-term damage. Due to the high rates which we were forced to use in order to reach 10-year equivalent LHC 
fluences, any determination of short term damage would not have been relevant for LHC conditions, anyway. 

Since no 7-irradiation data of recent CMS crystals existed at dose rates corresponding to those caused by the proton 
beam, we made comparative irradiations at an intense ^°Co facility. 

Only at proton fluences well beyond 10^^ cm^^, i.e. doses much higher than 1 kGy, we start to observe clear differ- 
ences in the LT-characteristics of crystals exposed to protons and those exposed to photons. The main differences 
are 

1 . The proton induced damage increases linearly with fluence while the 7-damage saturates after a few hours 
of exposure at 1 kGy/h. 

2. No correlation is seen between the ijij^^{420), from the radiation hardness pre-characterisation with ^°Co 
photons, and the proton induced /iiND(420). In the case of our 7-irradiated crystals this correlation is nearly 
perfect. 

3. In proton-irradiated crystals, the band-edge shifts towards longer wavelengths and /iiND(A) c>c A^^, while 
the band-edge remains stable in the 7-irradiated crystals and a clear absorption-band structure appears in 
A*ind(A). 

^' In the simulations for the EE, one crystal within a big matrix was considered, while for the beam test a single crystal in 
infinite void was simulated. 
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4. After ^ $p=5. 4x10^"^ cm^^, i.e. ^ 1 00 kGy, we observe ^ind (420) wl5m^^ and a linear rise with fluence. 
For a crystal with similar radiation-hardness characteristics, we observe after a 7-dose of ^50 kGy only 
Mind (420)=0.3 m^^ and no indication of further increase towards higher doses. 

It was not possible to raise the 7-intensity beyond 1 kGy/h, in order to obtain a comparison with our highest proton 
beam intensities. It is usually claimed that the saturation level of 7-damage increases with dose rate. For protons, 
however, we observe no difference between rates of 10^^ cm~^h^^ and 10^"^ cm^^h^^. 

A thorough exploration of fitting various functions to the recovery data showed, that Eq.|3lis suitable to fit all data 
accumulated so far In particular we discovered that the quality of the fits does not suffer, if the time constants, 
Ti, are universal with values of 17.2 days and 650 days. Only the absolute amount of damage is very different for 
the two radiation types, which is reflected in a different dose dependence for the fitted amplitudes Ai, as shown in 

Fig-ini 

Our results provide strong evidence that proton exposure of PbW04 crystals causes damage that cannot be repro- 
duced with 7-irradiation. This damage seems to be independent of rate and increases linearly with accumulated 
fluence. We strongly emphasise that the difference we observe starts to manifest itself only at proton fluences well 
beyond 10^^ cm^'^. Below this fluence the "normal" 7-like damage, due to minimum ionising tracks, is dominating 
to such an extent that it is probably impossible to experimentally see any difference between protons and photons. 

The observation, that /iiND(A) of proton damaged crystals is consistent with Rayleigh scattering from objects 
much smaller than 0.4 /im, supports our assumption that the specific proton damage is due to the very dense local 
ionisation of fragments created in inelastic nuclear collisions. In this case it is likely that the 20 Ge V/c proton beam 
used in our tests is somewhat too pessimistic with respect to the hadron spectrum expected in the CMS ECAL. 
Our simulations indicate that the damage ratio between 20 Ge V/c protons and the CMS ECAL hadron spectrum, 
per unit star density, might be about a factor of 4. But this remains to be verified by further irradiations with 
low-energy pions, which provide a better representation of the ECAL conditions. 

A Induced radioactivity 

The creation of nuclear fragments during hadron irradiation causes the crystals to become radioactive. Since 
PbW04 includes the heavy elements W and Pb, the 20 Ge V/c proton beam is able to produce a huge number of 
different nuclides. Most of these are either stable or short-lived, but a fraction of them have long half-life and 
are responsible for the remnant radioactivity of the irradiated crystals. According to simulations, there are about 
1800 different nuclide species with an activity above 1 Bq/cm^ at the end of an exposure to 10^^ p/cm^. Even after 
50 days of cooling, 150 nuclides remain above this activity threshold. 

The predictive simulation of crystal radioactivation proceeded in three steps. First, during the full FLUKA simula- 
tion described before, all residual nuclides produced in the crystal were recorded. This stage relies on the accuracy 
of the FLUKA models to describe the hadron-nucleus interaction in full detail. Earlier dedicated work has shown, 
that an accuracy better than a factor of 2 can be expected for heavy materials |30 311. In a second stage the evo- 
lution over time of the produced nuclide inventory was calculated with the DeTra code [32 1. This calculation is 
analytical and thus exact within the limits of accuracy given by the available nuclear data files. DeTra produces 
a table of all nuclide activities as a function of cooling time. For the final phase of the simulation these activities 
were converted into an emitted 7-spectrum which was then used to calculate Bind at our measurement point by 
integrating over the entire crystal and weighting with the simulated star density within the crystal. The integration 
was performed via a detailed photon transport with FLUKA. 

Figurel20l compares simulated -Djnd and all our measurements. For this purpose all measurements have been 
rescaled to ^ <i>p = IG^'^ cm^^, using the total proton fluence, as determined from activation of aluminium foils. 
The two curves in Fig.|20| corresponding to different beam intensities but the same integral ^ $p, show that by 
the time our measurements could begin only the value of ^ <I>p matters, not the rate. Thus a rescaling by ^ $p is 
justified for our comparison. 

In the small insert of Fig.|20lthe doses from individual crystals are compared with the calculated Dind- We observe 
a systematic underestimate in the simulation. On average the deviation amounts to 29% and even at worst the dis- 
agreement remains below a factor of 2, which still can be considered a very good agreement for such a prediction. 
We also have to recall two uncertainties in the measured values with respect to the simulation, namely the unknown 
efficiency of the Automess 6150AD6 above 1.3 Me V, which we probably underestimate, and a slight uncertainty 
in the size and position of the sensitive element within the device. 



20 




-5 -4 -3-2-1 2 3 4 

10 10 10 10 10 1 10 10 10 10 

Time (days) 



Figure 20: Comparison of simulated and measured radioactivity as a function of cooling time. All data have been 
renormalised to ^ "I>j,=10^'^ cm^^h^^ using the fluence information from Al foils, as given in the ^ $p column 
ofTable|3l 
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7.67 


10.4 


8.64 


Calculated 


5.09 


6.58 


5.48 


1.58xCalc/Meas 


1.05 


1.00 


1.00 



Table 5: Measured and calculated Di^d (/iSv/h) close to the front, in the middle and close to the rear of crystal a" 
(see Table|3}. The values correspond to 177 days cooling after the last irradiation. 

We used the -Dind measurements also to cross-check the longitudinal star density simulations shown in Fig.|3] 
For this purpose we measured Dind, in addition to the middle, also close to the front and back of crystal a", 
which had been exposed to ^ $p=5.4x 10^^ cm~^. The results are given in Table|5l Correcting for the systematic 
underestimate of the simulations in the case of a ", the relative front-middle-backratios are very well reproduced'^^ 
This good agreement increased our confidence that the simulations, summarised in Figs.|3]and|3 describe well the 
conditions within the crystal. 
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Figure 21: Correlation of measured Dmd of the crystals and $p as determined from ^^Na activity in the Al foils. 

B Technical aspects of $p determination 

The most reliable standard way to determine the fluence of the incident proton beam ($p) is by means of aluminium 
activation measuring the production of ^''Na. This is the isotope for which cross sections are known most reliably. 
However, the cross section Al(n,a)^''Na above ^6MeV is much larger''^' than the Al(p,X)^*Na cross section at 
20GeV (8.7 mb). Since the crystal produces a large flux of neutrons above this threshold energy, ^^Na is not a 
useful isotope to monitor $p, on top of the fact that its 15 h half-life makes it impractical for our purpose. A more 
suitable isotope is ^^Na. Due to its long half-life it integrates the fluence, and is insensitive to short-term beam 
fluctuations. 

The cross section for Al(p,X)^^Na at ~20GeV is not known as accurately as for Al(p,X)^^Na. Data in 1341 exhibit 
a significant spread but allow us to estimate 11.2±1.6mb. The ratio ^^Na/^^Na in a foil, irradiated in the absence 
of a crystal, gave 10.4 mb. Since this is based on the much better known production cross section of ^*Na, we use 
10.4 mb for the fluence analysis. However, we have to point out that even without a crystal, backscattered radiation 
in IRRADl increases the production of ^"^Na more than of ^^Na. Thus the 10.4 mb is likely to represent a lower 
bound. Consequently, "I>p is likely to be rather over- than underestimated. 

While certainly smaller than for Al(n,a)^*Na, the cross sections for Al(n,X)'^^Na are not known well enough to 
allow for an accurate analysis of the neutron contribution. With our best knowledge of the cross sections we obtain 
from the FLUKA simulation a 1 1 % contribution to the ^^Na activity from non-beam particles. This we propagate 
into our estimate of the incident proton fluence. 

The aluminium foils were cut to a size of 24 x 24 mm^, which roughly corresponds to the average crystal dimen- 
sions. Thus the foils were averaging the whole fluence incident on the crystal. In a single crystal, however, protons 
close to an edge contribute on average less star density, thus damage and activation, because a larger fraction of the 
hadronic cascade escapes, when compared to centrally impinging ones. In fact, the leakage of the cascade causes 
that the radiation intensity drops towards the edges even with a uniform beam, for which the star density in the 
centre is about 17% higher than the average over the whole crystal. However, would the crystal be surrounded by 

''^ The values in front and back are particularly sensitive to the exact longitudinal position and the assumed size of the sensitive 
element in the 6150AD6. 

The Al(n,a)^^Na cross section exceeds 100 mb between 10-30 MeV and then decreases towards higher energies 1331 . 
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other crystals, also covered by the beam, then these would compensate for the leakage and the profile would be flat 
and the total star density higher. Thus the relation between <l>p, determined with the aluminium foil, and the star 
density depends - in addition to the beam characteristics - also on the size and surroundings of the crystal. 

The activation of the crystal by inelastic nuclear interactions is roughly proportional to the star density inside the 
crystal. Unfortunately, this does not provide an absolute determination of <i>p, because star density and .Dind can 
be related to each other only through simulations, which introduce too large an uncertainty, as discussed in detail 
in appendixIXI However, since for all our irradiations the conditions within the crystals are the same, the effects 
of simulations cancel out if we consider only the relative differences between irradiations. Thus we expect a near 
to perfect correlation between -Di„d and the ^^Na activity of the aluminium foils. The only deviations could arise 
from significant non-uniformities of the beam spot or from a severe misalignment of the crystal during irradiation 
and, to a lesser extent, from the decay corrections, for which recourse to simulation is unavoidable. This correlation 
between $p and Bind, shown in Fig.|^ provided a very valuable cross-check of our fluence measurements. The 
Dind measurements were taken at different times, but all values have been corrected to 50 days after irradiation 
using the simulated decay of radioactivity in the crystals. As expected, the correlation is almost perfect, indicating 
that the activation foils give reliable and consistent fluence values. There are no signs of severe anomalies in any 
individual irradiation. 

Due to a faulty magnet in the PS, the early irradiations took place with a 20GeV/c proton beam while later 
exposures (crystals a" and h) were performed at 24GeV/c. The energy dependence of the Al(p,X)Na^^ cross 
section is negligible in this regime. However, as shown in Fig.|3] the 20% increase in proton momentum leads 
to more intense showering and thus to a higher radiation field within the crystal. The effect is very similar for all 
hadronic radiation quantities plotted in Fig.|3]but we use the 8.6% increase in star density, averaged over the crystal 
length. The LT measurements also average along the whole crystal and thus justify a simple correction by simply 
rescaling <i>p, as determined from the aluminium foils, by a factor of 1.086 for the two crystals irradiated with 
24GeV/c protons. In the case of ionising dose the effect is larger because the electromagnetic cascade is much 
better contained even in a single crystal. 



C Influence of beam non-uniformities 
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Figure 22: Transverse star-density profiles across a crystal for a perfectly uniform beam and for a beam with 5 
times lower intensity at x=-12 mm than at a;=12 mm. 

In order to study the effect of non-uniformities of the beam, we assumed that the beam intensity drops linearly 
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Figure 23: Effect of beam non-uniformity in case of crystal E. Compared are the I^ind at 4.5 cm distance from each 
long side and the /kind (420)-values in the centre and close to each side. The measured values are labeled with 'm' 
and the simulated ones with 's'. 

from one side of the crystal to the other'^). Crystal E was identified as having received the most non-uniform 
exposure. We adjusted the slope such that the D\nd values measured on E were roughly reproduced and ended 
up with a factor of 5 intensity difference across the crystal front face. A full FLUKA simulation was done with 
such a severe non-uniformity. Fig.|22l compares the lateral profiles of star density obtained from this simulation 
and the normal uniform one. Although a larger fraction of beam hits close to the edge, the total number of stars 
differs only by a few percent. Due to cascading, the variation in star density is much smaller than the factor of 5 
used as incident beam intensity. We also observe that in the centre of the crystal the effects of non-uniformity are 
at a minimum, which suggests that, provided the crystal alignment is the same in all LT measurements, even for 
the most severe beam non-uniformity encountered the influence on the results remains negligible in comparison to 
other uncertainties. 

The factor of 5 beam intensity difference used in the simulations of a non-uniform beam, as described before, was 
adjusted so that the relative^) variation of -Dind observed on the most non-uniform crystal, E, was reproduced when 
the simulated data were used to calculate Dmd on each side. 

The same simulated data were also used to integrate the star density along the light paths of the LT measurements. 
Assuming that stars cause the damage, the fit of Fig. 1151 allows to translate these star densities into /^ind(420) 
values, which can be compared with the measured values on absolute scale. 

Fi gurel23l shows such comparisons for I)i„d and /iiND(420) in the case of crystal E. The systematic underestimate 
of -Dind is consistent with Fig.|20| where the simulated £)i„d-values for crystal E are also slightly lower than the 
measurements. However, it must be emphasised that for such a simulation an agreement better than a factor of 
two, is remarkable. 

The simulated ^ind(420) values are in perfect agreement with measurements, except for a slight up-down asym- 
metry. This is most probably related to a beam non-uniformity in the vertical direction, while the simulated beam 
intensity was varied only horizontally. 

A determination of the beam profile using OSL (Optically Stimulated Luminescence) films showed that the real non- 
uniformities are not linear but exhibit a bump-like maximum on one side. 

The absolute value of ^ "I>p is always derived from the aluminium foils. 
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In relative terms, the non-uniformities in both /^ind(420) and Dind are very well reproduced by just considering 
the spatial variation of star density within the crystal. 
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